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This talk presents results of weak matrix elements calculated from simulations done on 170 32^ X 64 lattices 
at P — 6.0 using quenched Wilson fermions. We discuss the extraction of pseudoscalar decay constants fn, /k, 
fo, and /ds, the form-factors for the rare decay B —f K*y, and the matrix elements of the 4-fermion operators 
relevant to Bk, By, Bg. We present an analysis of the various sources of systematic errors, and show that these are 
now much larger than the statistical errors for each of these observables. Our main results are fn = 186(29) MeV, 
fo, = 224(16)Mel/, Ti = Ta = 0.24(1), BK{NDR,2GeV) = 0.67(9), and B8{NDR,2GeY) = 0.81(1). 



1. TECHNICAL DETAILS 

We briefly state the details common to all three 
quantities discussed in this talk and refer the 
reader to for details. Preliminary results 

based on 100 lattices were presented at LAT- 
TICE94 Q and the final analysis will be pre- 
sented elsewhere 

We calculate wall and Wuppertal source quark 
propagators at five values of quark mass given by 
K = 0.135 (C), 0.153 (S), 0.155 (C/i), 0.1558 {U2), 
and 0.1563 (C/3). These quarks correspond to 
pseudoscalar mesons of mass 2835, 983, 690, 545 
and 431 MeV respectively where we have used 
1/a — 2.33 GeV for the lattice scale. We con- 
struct three types of correlation functions, Wup- 
pertal smeared-local {Tgi^) and smeared-smeared 
(Fgs), and wall smeared-local {Twl)- The three 
Ui quarks allow us to extrapolate the data to 
the physical isospin symmetric light quark mass 
m = (m„ -|- md)/2, while the physical charm 
mass is taken to be C. The physical value of 
strange quark lies between S and Ui and we use 
these two points to interpolate to it. For brevity 
we will denote the six combinations of light 
quarks f/iC/i, U,U2, U1U3, U2U2, U2U3, UsU^hy 
{UiUj} and the three degenerate cases by {UiUi}. 
The FsL, Tss, and 3-point functions have been 
evaluated at the 5 lowest lattice momenta, i.e. 
p=(0,0,0),(l,0,0),(l,l,0),(l,l,l),(2,0,0). 



* Based on talks presented by Rajan Gupta and Tanmoy 
Bhattacharya. These calculations have been done on the 
CM5 at LANL as part of the DOE HPCC Grand Challenge 
program, and at NCSA under a Metacenter allocation. 



Renormalization Constants: We use the 

Lepage-Mackenzie tadpole subtraction prescrip- 
tion Its implementation consists of three 
parts in addition to writing the perturbative ex- 
pansions in terms of the improved coupling ay. 
One, the renormalization of the quark field y/Z^ 
changes from 



'2k \/8kc\/1 — 3k/Akc', second, 
the perturbative expression for 8kc in Z^p is com- 
bined with the coefRcient of ay in the one loop 
matching relations to remove the tadpole contri- 
bution, and finally the typical momentum scale 
once the tadpole diagrams are removed is taken 
to be g* = 1/a, i.e. both the scale at which a„ is 
evaluated and the scale at which lattice and con- 
tinuum theories are matched is set to q* — 1/a. 
We label this scheme TADl for brevity. The 
difference in results for q* = 1/a and q* — ir/a is 
used as an estimate of systematic errors due to 
fixing q*. Our data gives Kc = 0.157131(9) 0. 

Setting the quark masses: In |^ we show that 
a non-perturbative estimate of quark mass TO„p, 
calculated using the Ward identity, is linearly re- 
lated to (1/2k— 1/2kc) for light quarks, so either 
definition of the quark mass can be used for the 
extrapolation. We choose to use ninp, and fix m, 
nis and rric as follows. To get m we extrapolate 
the ratio M^/M^ to its physical value 0.03182. 
We determine nis by extrapolating M^/Mp to 
m and then interpolating in the strange quark 
to match the physical value. We find a ~ 20% 
difference between using M|./MJ or M^/Mp to 
fix nis, which we use as an estimate of the sys- 
tematic error. For we use k = 0.135 as we 
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have simulated only one heavy mass. With this 
choice the experimental values of AId, Md* and 
Md, lie in between the static mass Mi (mea- 
sured from the rate of exponential fall-off of the 
2-point function) and the kinetic mass defined as 
M2 = {d"^ E / dp\=o)~^ . The difference in final 
quantities between using Mi and M2 is taken to 
be an estimate of the systematic error in fixing 
rUc- 

The lattice scale a: To convert lattice results 
to physical units we use a{Mp). As discussed 
in iQ], the Mp data show a small but statis- 
tically significant negative curvature. We get 
l/a = 2.330(41) GeV from a linear fit to {C/^C/j} 
points, 2.365(48) GeV including a m'^/^ correc- 
tion term in the fit to the 10 UiUj, SUi, SS points, 
and 2.344(42) GeV including a term. Since 
all three estimates are consistent and the form 
of the chiral correction cannot be resolved we 
use the result from the linear extrapolation and 
assign 3% as an estimate of the systematic error. 

2. DECAY CONSTANTS 



The pseudoscalar decay constant fps is given 



by 



ZA{0\A'r^'\n{p) 



(1) 



where Za is the renormalization constant con- 
necting the lattice scheme to continuum A/5. 
We study, in addition to the 2-point correlation 
functions T. two kinds of ratios of correlators: 
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Using either tt or ^4 for the smeared source J 
gives 4 ways of extracting fps- Two more ways 
are gotten by combining the mass and amplitude 
of the 2-point correlation functions, i.e. {A/^P)ls 
and {PP)ss, and {AiAi)LS and (^4^4)55- 

The data satisfy the following consistency 
checks: the six ways of calculating fps described 
above, and at each of the five values of momen- 
tum, give results consistent to within 2cr . (The 
one exception is the p = (2,0,0) case where the 
signal is not good enough to ascertain that we 
have fit to the lowest state.) Even though these 



Figure 1. Plot of the Bernard-Golterman ratio R 
for the quenched theory. 
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estimates are correlated, consistent results do in- 
dicate that fits have been made to the lowest state 
and reassure us of the statistical quality of the 
data. We use the p = (0, 0, 0) data in our final 
analysis as it has the best signal. 
Quenched approximation When analyzed in 
terms of chiral perturbation theory (CPT), there 
are two consequences of using the quenched 
approximation. One, the coefhcients in the 
quenched theory are different from those in full 
QCD and uncalculable, and second, Sharpe and 
collaborators and Bernard and Golterman § 
have pointed out that there exist extra chiral logs 
due to the 77' as it is also a Goldstone boson in the 
quenched approximation. These make the chi- 
ral limit of quenched quantities sick. To analyze 
the effects of 77' loops Bernard and Golterman Q 
have constructed the ratio R = /i2//ii'/22' ap- 
plicable in a 4-flavor theory where mi = toi' and 
TO2 = . The advantages of this ratio in com- 
paring full and quenched theories is that it is free 
of ambiguities due to the cutoff A in loop inte- 
grals and O(p^) terms in the chiral Lagrangian. 
CPT predicts that 
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Figure 2. Plot of R for the full theory. 
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where 6 = m§/247r^/^ parameterizes the effects 
of the 7]', and Xguenched and Xjuu are given in 
||. At LATTICE94 the preferred fit (with 100 
configurations and no (mi— 7712)^ term) was to the 
quenched expression which gave S = 0.10(3) 0. 
The need for including the (mi — TO2)^ correction 
is shown in Figs. and ||. The fit to the quenched 
expression gives S = 0.14(4), however, based on 
X^, the fit to the full QCD expression is preferred. 
The caveat is that the intercept is 1.69(45) rather 
than unity. Thus, we cannot resolve the effects of 
77' from normal higher order terms in the chiral 
expansion, and neglect both in our analysis. 

Extrapolation to the physical quark 
masses: The data, shown in Fig. ^, indicates 
a break in the vicinity of m^ between the non- 
degenerate {SUi} and degenerate UiUi mesons at 
the la level, but no such break between the UiUi 
and the U2U2 cases. We thus use {UiUj} points 
to extrapolate to /^. Note that even though the 
slopes for the two fits to {UiUj} and {SS,SUi} 
combinations are different, the values after ex- 
trapolation are virtually indistinguishable. In 
Fig. ^ we show the extrapolation for heavy-light 
mesons for three cases (C, S, Ui) of "heavy" 
quarks. In all three cases we use a linear fit to 
the three Ui points for extrapolation to rn as de- 
viations from linearity are apparent if the "light" 



Figure 3. Plot of data for /^r versus m„p. The 
linear fit (solid line) is to the six {UiUj} points, 
with errors shown by the dotted lines. The dash- 
dot fine is a hnear fit to the four {SS, SUi} points. 
The vertical line at m„p ~ represents m and the 
band at m„p « 0.04 denotes the range of mg. 
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Figure 4. Extrapolation of heavy-light fps to m 
for three cases, C, S, Ui, of "heavy" quarks. The 
linear fits are to the three "light" Ui quarks, and 
the fourth point (light quark is S) is included to 
show the breakdown of the linear approximation. 
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quark mass is taken to be S as shown by the 
fourth point at r7i„p = 0.076. To get fx we in- 
terpolate to nis the result of the extrapolations of 
SUi and UiUi points to to. For fjj we extrapolate 
the three CUi, and for fjj^ we simply interpolate 
between CUi and CS points. 
Results at f3 — 6.0: Our final results using 
TAD! scheme along with estimates of statistical 
and various systematic errors are given in Table |l|. 
From the data it is clear that systematic errors 
due to setting TOc, the lattice scale, and Za are 
now the dominant sources of errors. 
Continuum Limit: To extract results valid in 
the continuum limit we include data from the 
GFll iP = 5.7,5.93,6.17) [|lO|, JLQCD {f3 = 
6.1,6.3) [H, and APE (/? = 6.0,6.2) Q Col- 
laborations. We have attempted to correct for 
as many systematic differences, however some, 
like differences in lattice volumes, range of quark 
masses analyzed, and fitting techniques, remain. 

Assuming that lattice spacing errors are 0{a), 
a linear fit versus AIpU gives 

U/Mp = 0.156(7) (expt. 0.170), 
fx/Mp = 0.171(6) (expt. 0.208). 

with x^/dof = 1-6 and 1.7 respectively. The 
change from the GFll results is marginal as the 
fit is still strongly influenced by the point at 
f3 — 5.7, which may lie outside the domain of 
validity of the linear extrapolation. A linear ex- 
trapolation excluding the f3 = 5.7 data gives 

U/Mp = 0.170(14), fK/Mp = 0.187(11), 

with x^/dof = 2.1 and 1.9 respectively. Using 
ms{M^) would increase Jk by ~ 2%. Given this 
difference in the extrapolated value depending on 
whether the data at /3 = 5.7 is included or not 
makes it clear that more data are required to 
make a reliable a — > extrapolation. 

The fo and fo, data at /? > 6.0, in TADl 
scheme, and using ms{MK) are shown in Fig. ||. 
The APE collaboration use Afi for the meson 
mass. For consistency we have shifted their data 
to M2 using our estimates given in Table |^. A 
linear extrapolation to a = then gives 

fo = 186(29) MeV, /d, = 218(15) MeV, 



Figure 5. Extrapolation to the continuum limit 
of Jd and (in MeV) data. Our data is shown 
with the symbol octagon, the plus points are from 
the JLQCD Collaboration |1^, a nd the diamonds 
label the APE collaboration l|ll| data. 
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with I dof — 2.2 and 2.0 respectively. Using 
TOs(M0) increases fa, to 224(16) MeV. The qual- 
ity of the fits are, however, not very satisfactory. 
The bottom line is that in order to improve the 
estimates the various systematic errors that have 
not been included in the a —> Q extrapolations 
presented above need to be reduced. 

3. THE RARE DECAY B K*-f. 

We discuss the applicability of heavy quark ef- 
fective theory (HQET) and pole dominance hy- 
pothesis (PDH) to extract the form-factors Ti 
and T2 at = and mheavy = 'rri},. The 
technical setup is the same as described in |^ 
for the calculation of semi-leptonic form-factors, 
and the quality of the signal is similar to that for 
D ~* K*lv decays. 

PDH: states that the behavior of all form- 
factors is 

fm = f{Q)/{l-QVM^), (3) 

where M is the mass of the nearest resonance 
with the right quantum numbers. To test PDH 
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Table 1. Our final results at f3 ~ 6.0. All dimensionful numbers are given in MeV with the scale set by 
Mp. For the systematic errors due to mg, rric, q* we also give the sign of the effect. We cannot estimate 
the uncertainty due to using the quenched approximation, or for entries marked with a ?. 
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Figure 6. Three types of fits to test behavior 
of Ti using CU3 U1U3 transition. HQET sug- 
gests a dipole fit if one assumes PDH for T2. The 
data prefer a pole fit but with the resonance mass 
smaller than the lattice measured value Mn* used 
to plot the data. 
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we make two kinds of fits: (i) single parameter 
"pole" fit where M is the lattice measured value 
of the resonance mass, (ii) two parameter "best" 
fit where M and /(O) are free parameters. Typ- 
ical examples of these fits are shown in Figs. |^ 
and 0. Overall, Ti is well described by the "pole" 
form, whereas T2 has a "flat" dependence. We 
take the "best" fit values for our final estimates. 
HQET: To leading order in and in the mass 
of the heavy quarks, HQET implies (for heavy to 
heavy transitions) that the combinations 



ruB 



Ti(g2 



•■K 



m,B 



'K 



1 - 



-(4) 



are independent of the masses of the heavy quarks 
for fixed velocity transfer. Since Ti{Q^ = 0) = 
T2(Q2 = 0) for all m,, this HQET relation and 
the PDH, Eq. ^, cannot hold simultaneously. In 
fact, for heavy quarks (ms + niK*) ~ rripoie, 
therefore, if T2 fits the pole form then Ti must 
be a 'dipole'. Instead our data, as exemplified in 
Figs. H and |^, prefer a flat T2 and a pole behavior 
for Ti. 

Dependence on quark mass: Figures || and 
^ show examples of the variation of Ti(0) and 
T2(0) with quark masses. There is signiflcant de- 
pendence on the mass of the quark C decays into 
(which is a kinematic effect), and a slight depen- 
dence on nispectator resulting in the small increase 
in slope between CUi UiUi and CUi SUi 
cases, which is consistent with HQET. 
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Figure 7. Fits to test behavior of T2. 
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Figure 8. Extrapolation of T\{Q'^ = 0) to m„. 
The interpolation to for T\{B — + A'*7) is 
done using the points labeled by squares (s — S) 
and octagons (s = Ui), while Ti{B p^) is ob- 
tained by extrapolating the degenerate qq points 
(crosses). 
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Figure 9. Extrapolation of T2(Q^ — 0) to m„. 
The interpolation to rus for T2{B K*j) is 
done using the points labeled by squares (s = S) 
and octagons (s = Ui), while T2{B — > p7) is ob- 
tained by extrapolating the degenerate qq points 
(crosses). 
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Table 2. Estimates of form factors in 3 commonly 
used renormalization schemes defined in Q . The 
data satisfy Ti{Q^ = 0) = T2(Q2 = 0). The 
last four rows give T{Q^ = 0) extrapolated to to;, 
using the 4 methods discussed in the text. 
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0.39(2) 
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Table 3. Estimates in TADl scheme at = 0. 
The variations give estimates of systematic errors. 
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Extrapolation in nxheavy ■ The need to extrap- 
olate the results obtained at rriheavy ~ ttIc to m?, 
using HQET, Eq. ^, introduces a very large un- 
certainty as shown by the four ways of analyzing 
the data. Methods 1 and 2: we take the value of 
T2 at zero recoil extrapolated to m^ and m and 
scale it to m?, using HQET. We can then esti- 
mate the value at = assuming pole domi- 
nance holds for T2 at mf, (advocated by A. Soni 
at this conference), or by using a "flat" behavior 



as shown by data at m/jg 



Method 3 (4): 



Scale T2{Q^ = 0) (Ti) assuming the joint validity 
of HQET and pole dominance. This implies the 
scaling relations T2{Q^ = 0)m^^^^^ = constant 

and Ti((5^ = 0)m^g^^^ = constant. The results 
along with their variation with the tadpole sub- 
traction prescription, type of fit, mg, and the def- 
inition of heavy-light meson mass (Mi or M2) are 
shown in tables |^ and ^. 

Results at /3 = 6.0: Methods 1,2 and 3,4 reflect 
the same contradiction. The value is either 0.08 — 
0.1 or 0.23 — 0.25 depending on what we assume 
for the scaling behavior. With present data we 
assume that the flat behavior for T2 and pole 
dominance for Ti persists all the way upto the 
physical value of mf,. Then, using the best fit, 
TADl subtraction prescription, ms(M0), and Mi 
for the meson mass, we get Ti — T2 — 0.24(1). 
Further progress requires clarification of the 
behaviour of the form factors and an estimate of 



the violations of leading order HQET predictions. 

4. B-parameters 

We present an update on results for Bk^ 
Bj, Bg with Wilson fermions evaluated in the 
NDR scheme with TADl subtraction prescrip- 
tion. Note that both q* and the matching scale 
between the lattice and continuum theories are 
taken to be 1/a. Thereafter, the results are run 
to 2 GeV using the 2-loop relations, however the 
change is minimal. 

To analyze the lattice data (illustrated in Ta- 
ble ^) we consider the general form, ignoring 
chiral logs, of the chiral expansion of the AS* = 2 
4-fermion matrix elements with Wilson fermions 



5im\ 



+ S2mj^piPf + 5:i{piPf f + . . (5) 



This follows from Lorcntz symmetry as m^ and 
Pi ■ Pf are the only invariants. 
Bk' The terms proportional to a,/? and are 
pure lattice artifacts due to mixing of the AS* — 2 
4-fermion operator with wrong chirality opera- 
tors. To isolate these terms we fit the data for the 
lightest 10 mass combinations and for the 5 val- 
ues of momentum transfer using Eq. ^ as shown in 
Fig. |l^. (Similar values for the six coefficients are 
obtained from fits to the 6 lightest combinations.) 
We find that the three 6i are not well determined; 
only 62 is significantly different from zero. More 
important, the coefficients "(,62,63 contain arti- 
facts in addition to the desired physical pieces 
which we cannot resolve by this method. We 
simply assume that the 1-loop improved opera- 
tor does a sufficiently good job of removing these 
residual artifacts. The result then is 

Bk{NDR, 1/a) = -i+{62+ 63)Ml = 0.65(10). 

A second way of extracting Bk using Eq. || is to 
combine pairs of points at different momentum 
transfer: 

(EiBKiqi) - E2BK{q2))/{Ei - E2) = 

7 -I- 62m'^ + 5'im{Ei + E2). 

This procedure directly removes a, (3 and 5i but 
requires a correction to the 63m{Ei + E2) piece. 
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Figure 10. Six parameter fit to the Bk data. 
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for which we use the value of 6z extracted from 
the fit. The results of this analysis for the 10 
light mass combinations are given in the third 
column of Table ^. Interpolating to m.K , we get 
BK{NDR,l/a) = 0.66(9). 

Bk ■ The 2-loop running of Bk defines the renor- 
malization group invariant quantity Bk fl^ 



Bk = asifJ.)' 



-2//3o , 



1 



,{pi)J/47T)BK{fi) 



where J = (/3i7o - /3o7i)/2/3^ = 2.004 for Uf = 
0. Under running, Bk increases as ^ is de- 
creased. Thus BK{NDR,l/a) = 0.66(9) be- 
comes BKiNDR,2GeV) = 0.67(9), and Bk is 
0.90(14). For comparison, the Staggered value 
calculated at /3 = 6.0 is B^f (iVDi?, 2 GeV) = 
0.67 — 0.71 depending on the lattice operators 
used 14 1 . An update on staggered results and 
issues of extrapolation to a — > have been pre- 
sented by JLQCD at this conference p^ . 
Bd- The CS, CUi data show no significant 
variation with momentum transfer as shown in 
Table ^. The theoretical analysis of artifacts 
in heavy-light mesons has not yet been com- 
pleted; indications are that all 6 terms con- 
tribute. Therefore we simply extrapolate the CUi 
data to m to get BD[NDR,l/a) = 0.78(1) or 
Bd{NDR,2 GeY) = 0.79(1). 



Table 4. Bk, B-; and Bg in NDR renormalization 
scheme at matching scale — 1. 
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Br 



B. 
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CC 
CS 
CUi 

CU2 

CU3 

SS 

SUi 

SU2 

SU3 

UiUi 

U1U2 

U1U3 

U2U2 

U2U3 

U3U3 



0.92(1) 
0.83(1) 
0.81(1) 
0.80(1) 
0.79(1) 
0.56(0) 
0.44(0) 
0.38(1) 
0.34(1) 
0.24(0) 
0.11(1) 
0.00(1) 
-0.09(1) 
-0.29(1) 
-0.60(2) 



10.69(30) 
10.69(23) 
10.68(21) 
10.67(19) 
10.68(16) 
10.68(14) 
10.66(12) 
10.67(11) 
10.65(10) 
10.63(09) 



0.86(1) 
0.82(1) 
0.81(1) 
0.81(1) 
0.80(1) 
0.72(1) 
0.70(0) 
0.68(0) 
0.67(1) 
0.66(0) 
0.64(0) 
0.63(0) 
0.62(0) 
0.60(0) 
0.58(1) 



B7 and Bg: The chiral expansion is simi- 
lar to Eq. H with all 6 coefficients contain- 
ing artifacts and physical pieces. Ignoring 
the artifacts we get Br{NDR,l/a) = 0.60(1) 
and BsiNDR,l/a) = 0.81(1) or with 2- 
loop running Br{N DR,2 GeY) = 0.59(1) and 
Bs{NDR,2 GeY) = 0.81(1). (We find that Bg is 
very insensitive to changes in fi as the running of 
the matrix element is almost completely canceled 
by that of its vacuum saturation approximation, 
while in the case of Bj the two add.) Our es- 
timate of Bs is smaller than that used in the 
Standard Model analysis of e'/e |l^]. Since a 
smaller Bg means larger e'/e, the calculation of 
Bg, is important phenomenologically. Work is in 
progress to understand and remove various lattice 
artifacts and make our estimate more reliable. 
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